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Abstract: We study mesons in N = 4 super Yang-Mills theory with fundamental flavors 
added at large 't Hooft coupling using the gauge/gravity correspondence. High-spin mesons 
are well described by using semiclassical string configurations. We determine the meson 
spectrum at finite temperature and in a background magnetic field. 
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1. Introduction 

The determination of the bound state spectrum in a Coulomb potential is one of the oldest and 
most important successes of quantum mechanics. Since the fine-structure constant in nature 
is small the binding energy is small, E^ind ~ a 2 m << m. The bound states are essentially 
non-relativistic. It is interesting to understand what the properties of bound states look like 
in the opposite limit of extremely strong Coulomb interactions. One system that realizes 
such strongly coupled Coulombic bound states is M = 4 supersymmetric Yang-Mills theory 
with J\f = 2 supersymmetric flavors in the limit of large number of colors N c and large 't 
Hooft coupling. The N = 4 theory is conformal, which guarantees that the potential between 
external test particles is of the form V{r) = ——. In the large N c , large 't Hooft coupling 
limit the theory can be solved using the AdS/CFT correspondence [1-3]. In particular the 
coefficient c\ in the Coulomb potential between external test sources has been determined 1 
to be 

4vr 2 r- 

« - ^ (") 

from a Wilson line computation [4,5]. 

To realize Coulombic bound states, one needs to add Nf flavors of dynamical fundamental 
matter (quarks) to the M = 4 theory. In the large A^ c , finite Nf limit this can be achieved 
via adding flavor probe branes on the supergravity side [6]. In this case conformal invariance 

1 We are using conventions where A = 2<?y M iV = 4-ivg s N. 
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Figure 1: Quark/antiquark potential V(r) in units where the quark mass is m = -j^IO. 

alone doesn't completely fix the form of the potential. The quantum corrected Coulomb 
potential instead has the form 

V(r) = ——f(mr). 

where m denotes the mass of the fundamental flavors. At strong coupling, the function / can 
be determined from a static semiclassical string calculation in the bulk. This was first done 
in [7] . The resulting potential is depicted in figure (Q) . At large separation (alternatively, at 
large mass) we need to recover the result for infinitely heavy test quarks and hence / — * 1. 
At small separation, / ~ m 2 r 2 so that the potential becomes linear in distance, V(r) = T e ffr 
with 

Teff= 7T (L2) 

The system is highly relativistic in this regime: the glue forms flux tubes that are well 
described by a relativistic string. 

At parametrically small separations, mr ~ the determination of / in terms of a 
semiclassical string breaks down. To ensure that the mesons we study are insensitive to the 
details of the potential in this ultra-short distance regime we study mesons with large angular 
momentum, J » 1. For J » vA the quark/antiquark separation is in the Coulombic 
regime and the corresponding mesons reproduce hydrogen-like behavior. For 1 << J << y/\ 
the mesons are dominated by the linear potential and one can compare the spectrum to that 
of open relativistic strings in flat space. In all cases the meson masses are of order the quark 
mass. On top of this spectrum of highly spinning strings, there are deeply bound mesons with 
J = 0, 1 that are dual to the small fluctuations of the flavor brane worldvolume fields [7] . Their 
masses are of order and they have been analyzed before at finite temperature [8-12] and 
non- vanishing magnetic field [13-15] or finite chemical potential [16]. For a review see [17]. 
Our studies complement these results for J = 0, 1 supergravity mesons by analyzing the 
highly spinning mesons with semiclassical strings in AdS. 

In the next section we start with a review of the meson properties at zero temperature. 
In section 3 we then analyze mesons and their dissociation at finite temperature. The results 
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are very similar to those obtained in [18] for the closely related Sakai-Sugimoto model [19]. 
In section 4 we study mesons in a background magnetic field and compare to results for 
the Zeeman splitting both in hydrogen and for a relativistic string in the appropriate limits. 
Finally, in section 5, we bring the two together and study the dissociation of mesons at finite 
temperature in the presence of a magnetic field. 



2. Meson properties at zero temperature 

Highly spinning mesons at zero temperature were first studied in [7]. The corresponding 
bulk configuration is a string with both of its ends on the flavor brane. For large angular 
momentum J >> 1, the length of the string is much larger than l s and it can be analyzed 
using the classical Nambu-Goto action in an AdS background. Writing the background AdS 
metric (setting the curvature radius to 1) as 

dr 2 

ds 2 = -h(r)dt 2 + r 2 dx 2 + -f-- (2.1) 

h(r) 

with h(r) = r 2 we are interested in stationary string configurations whose endpoints describe 
a circular orbit in the x directions. Writing dx 2 = dp 2 + p 2 d9 2 we are hence looking for a 
string whose worldvolume is described by 

9 = cot, r = r(a), p = p{a). (2.2) 

The Nambu-Goto action for this ansatz becomes 



S=-^- J dtda^(r' 2 + p' 2 r 2 h(r)) (l - = J dtda^J (1 — uj 2 p 2 ){r' 2 + r A p' 2 ) 

(2.3) 

where we used that in the units where the AdS curvature radius is one, the string tension is 
simply 

A = ^ (2.4) 
2jra' 2vr v ; 

The energy and angular momentum of the string follow immediately 



8L f 2 r> 2 + r*p> 2 

= ~8uj = ~2tt J ^ y ~Y^uj 2 ~p 2 ~ ^ 

dL VX f, l r> 2 + r*p> 2 

The equation of motion is straightforward. The Neumann boundary condition on p that 
follows from the variation of this Lagrangian is 



p Op' 2vr up V r' 2 + rV 2 



- 3 - 



which simply demands p' = at the end, that is the string ends on the flavor brane at a right 
angle at both ends. Since the configurations we are interested in are even under reflection 
around the string midpoint, one can alternatively specify a boundary condition each at the 
string midpoint and at a single endpoint. 

To proceed, one can further fix the a reparametrization invariance. Two convenient gauge 
choices seem to be either the a = p or the a = r static gauge. Both lead to large gradients 
which makes the numerics unstable. For the numerics we found that a very efficient choice is 

a = p + r (2.8) 

and hence p' +r' = 1. Numerically, we impose boundary conditions in the IR, demanding that 
a turnaround point exists at p = 0. That is, at a = ro we impose the boundary conditions 

r(r ) = r , r'(r ) = =>- p(r ) = 0, p'(r ) = 1. (2.9) 

For every choice of turnaround position ro one then looks for the first r m > ro at which 
the UV boundary condition p'(r m ) = is satisfied. This way one avoids shooting and every 
numerical run will generate a valid solution for some mass m that is given in terms of r m by 

m = ^—r m . (2.10) 

With this numerical procedure we were able to reproduce the spectrum of mesons E{J) that 
was previously obtained in [7]. Our numerical results, together with the analytic limits we are 
about to discuss, are displayed in figure @. We use the same gauge and numerical methods 
for our analysis in later sections. 

As discussed in the introduction there are two regimes in which we expect to be able 
to compare the meson spectrum against analytic formulas. For large J » vA we are in a 
Coulombic regime where the potential V(r) = —c\jr with c\ given in eq.( |l.l| ). The system 
is non-relativistic in this regime. The standard non-relativistic boundstate spectrum in a 
Coulomb potential is given by (including the rest-mass 2m of the two quarks) 

, . . c? m* 47r 4 A 

E(J)/m = 2-^- = 2- fw - 2 (2.11) 

where we used that the reduced mass is m* = ™^]^ 2 = -y for two equal mass particles. The 
other tractable region is 1 << J « y/\ in which we expect standard Regge behavior 

E 2 = 2irr e ffJ (2.12) 



where the effective tension is given by eq.(1.2) so that 



E/m = 2vr^J. (2.13) 
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E/m 




Figure 2: Meson energy E/m as a function of x/j/A 1 / 4 . E/m = 2 is the rest mass of the quarks, 
so bound states only exist for E/m < 2. Also indicated are the two limiting case of a pure Coulomb 
potential which is a good approximation for J >> VA and Regge behavior for 1 << J << Va. 



3. Meson dissociation at finite temperature 

In this section we study how the spectrum of highly spinning mesons is modified at finite 
temperature. In particular, we study the dissociation of the most weakly bound mesons 
as the temperature is raised. The disappearance of these bound states leads to a critical 
J cr it(m/T), where mesons with J > J cr n become unstable. The mesons of flavored N = 4 
have been proposed as a model for charmonium boundstates at finite temperature e.g. in 
reference [20] where their dissociation as a function of linear velocity has been analyzed. The 
study of meson dissociation is then potentially relevant for studies at heavy ion colliders. 
Meson dissociation in the closely related Sakai-Sugimoto model has been analyzed in [18] 
with qualitatively similar results. 

In the numerical analysis, the only difference with the zero temperature case is that in 

4 

the action we need to set h(r) = r 2 — with r/j = ttT instead of the h{r) = r 2 we had before. 
In interpreting the result, some care has to be taken in what one means by mass. As discussed 
in [21] at finite temperature (and the same comment will later apply to finite magnetic field) 
there are several notions of mass which all coincided in the B = T = case. One is the mass 
parameter that appears in the Lagrangian, m^. This is related to the position r m at which 
the flavor brane ends in a non-trivial fashion as first discussed in [8]. One needs to know the 
exact embedding of the brane, which changes both in response to the temperature and the 
magnetic field. These embeddings have been obtained in [15,22] A physical definition of mass 
that is directly related to r m is the rest mass of a quasi-particle, which was denoted M res t 
in [21]. This is simply the energy that one has to pay to create a flavored quasi-particle and 
hence just given by the length of the string stretching from the horizon to the flavor brane 
times the tension. For our purposes, we will continue referring to m = as "the mass". 

In the language of [21] this is M res t + Am, latter being the thermal mass shift, Am = . 
All other notions of mass can simply be determined in terms of m, for details see [21]. With 
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Figure 3: J cr u as a function of ^. 

this definition of mass the minimal mass is M res t = and hence m = Am = This is 

the point where the brane touches the horizon. This minimal value of m appears in several 
of our plots. 

Figure @ displays our result for the critical angular momentun as a function of m/T. As 
in the zero temperature case, one might be tempted to give an interpretation of J cr it(m/T) 
using a field theory model of two massive quarks moving in an effective potential. However, 
such an interpretation is only sensible in the non-relativistic regime. At finite temperature, 
due to screening, the quark/anti-quark potential becomes short range and the non-relativistic 
regime gets removed first. For relativistic bound-states, a significant fraction of the angular 
momentum and energy gets carried by the flux tube for which we only have a 5d descrip- 
tion. Qualitatively, the fact that J cr it drops as the temperature rises is consistent with an 
interpretation of bound states melting. 

Last but not least note that our semiclassical string configurations only probe the short 
distance part of the quark/anti-quark potential that scales as y/X. As argued in [23] there 
will be an exponential tail on top of this with an order one coefficient. All we can say with 
confidence is that for angular momenta above J cr it mesons with binding energies of order y/X 
do not exist. There might be bound states with order one binding energies which would be 
represented by spinning strings stabilized by quantum corrections. 

4. Zeeman effect 

4.1 General Considerations 

In this section we study the effect of a background magnetic field on the meson spectrum. At 
zero magnetic field rotational symmetry ensures that the energy of the mesons only depends 
on the magnitude J of the angular momentum and not its direction. As usual, the magnetic 
field breaks this degeneracy and introduces additional dependence of the meson masses on the 
scalar product J-B = rrijB. With this definition of rrij, it takes values between — J and +J as 
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usual. Since we are working in the limit of large J we do not see the quantization of rrij. For a 
generic orientation of J and B, the orbits of the quark/antiquark pair will no longer be circular 
and our ansatz for the bulk string configuration is too simple-minded. However, the same 



ansatz eq.(2.2) can still capture the two extreme cases of rrij = ±J, that is the magnetic field 
is either completely aligned or anti-aligned with the angular momentum. These two extreme 
cases are sufficient to analyze the splitting in the magnetic field background. For small B, 
one can linearize the energies in J ■ B and find the standard behavior 

AE = a (J) rrijB. (4.1) 

Calculating the energy splitting for rrij = ±J is completely sufficient to determine a(J). 
When non-linearities in B become important, we can still exactly trace the two extreme 
levels. Values of — J < rrij < J are expected to lie in between. They need no longer be 
equally spaced. 

There are actually two cases we will discuss and need to distinguish: quark/anti-quark 
pairs with equal opposite charges or with equal charges. Of course they will always be 
oppositely charged under the gauged SU(N) symmetry. For a single flavor brane, there is 
just one global U(l) gauge field living on the flavor brane, baryon number. The quark and 
the anti-quark carry equal and opposite charge under this global symmetry. Since quark and 
anti-quark have the same mass, the orbital magnetic moment of the two actually cancels and, 
like in positronium, the linear Zeeman splitting does not affect states with different values of 
rrij. To see any non-trivial linear Zeeman splitting we hence also study the case of equally 
charged string endpoints. One way to achieve this is to look at 2 flavor branes A and B, so 
that the global flavor symmetry is U(1)a x U(1)b an d then turn on equal opposite fields in 
the two U(l) factors. While the AA and BB strings will, like before, have equal opposite 
charges under this field, the AB or BA strings will have equal charges. For simplicity, we only 
look at 2 flavors of equal mass, but obviously our method can be extended to the case of two 
unequal mass flavors. The corresponding field theory and gravity configurations are sketched 
and contrasted in figure (H). 

In the case of equal opposite charges, there is no linear Zeeman splitting since the mag- 
netic moments of the two particles cancel. There can however be a correction to the energies 
quadratic in B. Since the only two configurations we can realize in the supergravity corre- 
spond to J ■ B = ±JB, in this case there will be no splitting and the two configurations will 
have identical energies. It is indeed easy to see both in the field theory and the supergravity 
that in the case of equal opposite charges, there is a symmetry that takes J — ► — J and 
at the same time exchanges the endpoints of the string. This symmetry operation takes one 
configuration into a different one with identical energy but opposite sign of J -B. For the case 
of 1 << J « vA we can compare our supergravity analysis with the study of a Nambu-Goto 
string in a background magnetic field. This was analyzed in the classic paper [24] where it 
was found that in the case of equal opposite charges there is no change in the string spectrum 
whatsoever, even though the actual string configuration changes: the string starts to fold 
back onto itself as indicated in the left of figure (|). So for J << a/A the energy shift 
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Equal Opposite No magnetic field Equal Charges 

Charges 

Figure 4: Sketch of supergravity and field theory configurations for the cases of equal opposite and 
equal charges of the string endpoints. In the field theory picture the circles indicate the position of 
quark and anti-quark, whereas the thick lines are the flux-tubes. In the case of equal opposite charges 
there is no net magnetic moment since the contributions from quark and anti-quark cancel. In this 
case there is a second configuration with the same energy that is the mirror image of the one depicted. 
In the case of equal charges the two contributions to the orbital angular momentum add up. There 
are two configurations with energies E± = E ± AE. 



function of B 2 has to vanish. For J >> \/A one expects a non-trivial shift in the energy, even 
though there will not be any splitting between the two configurations with opposite J ■ B. 

In the case of equal charges, we should see the familiar pattern of Zeeman splitting with 
a shift proportional to J • B so that for any given J and B we will find two configurations 



with circular orbit and energies E(B = 0) ± AE with AE of the form indicated in eq. (4.1). 
At J » a/A we can compare our results to an analysis of a purely Coulombic Hydrogen-like 
system. In this case we expect 

AE_„ B JB_ B j 



m m 



2m 2 



where Ms = 2m is the standard Bohr magneton for particles of unit charge. ^ and ^ are the 
dimensionless quantities measuring the energy and the magnetic field in units of the quark 
mass. Since our theory is conformal, the physics can only depend on these ratios and not on 
E, B and m individually. For 1 << J << \/A one can once more compare to the spectrum 
of a Nambu-Goto string with tension r e // in a background magnetic field. For equal charges 
it was found in [24] that for large J one has 

E 2 = 2vrr e// (l-e)J (4.3) 
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where 



2 B 2 B 

- arctan( ) = + 0(B 3 ) 

7T T e ff TTT eff 



so that to linear order in B we expect for the change in energy 



AE 



1 



m y/27TT e ff 

4.2 Supergravity analysis 



V7* 



B 



(4.4) 



(4.5) 



The presence of a constant magnetic field in the (p, 9) plane with F = Bdx Ady = Bpdr A d9 
only modifies the string action via the additional boundary terms. 



AS 



drA 



drA 



(4.6) 



where refers to the right and left boundary (which in our gauge choice is not simply at 
a = 7r, 0) and the upper (lower) sign corresponds to the case of equal opposite (equal) charges. 
Choosing A p = gauge we have as the only non-vanishing component 



A 6 = B P - 
2 



and the action ( |4.6| ) becomes (noting that 4f = oj) 



AS = B / dr 



2 

p OJ 



=F B dr 



p 2 oj 



(4.7) 



(4- 



This leads to two basic changes when compared to the analysis of section 2. First, the 
boundary condition gets modified to 



7T„ 



Fp8 @ 



-Bpoj. 



Second, the boundary terms in the action give an additional contribution to J = S-j 



AJ= B r — 
2 



=F B— 
2 



(4.9) 



(4.10) 



The corresponding contribution to the energy E = oj§^ — L cancels since the boundary action 
is linear in oj. A similar discussion appears in [25] where mesons in a non-commuative theory 
where analyzed which corresponds to a gravity solution with non-trivial background H-flux. 

For the case of two equal charges, the configurations are once again symmetric under 
reflection around the midpoint and one can then simply perform the numerical analysis along 
the same lines as in the zero temperature, zero B-field case. For every choice of turnaround 
point tq one generates a string configuration that solves the equation of motion. For every 
choice of r m one then can subsequently read off a B that allows one to have the boundary 
condition satisfied at that value of r m . The change in angular momentum in this case is given 



- 9 - 




Figure 5: Numerical results for E/m as a function ofB/m 2 for J/VA = 0.00159 and for J/y/X = 3.34 
together with analytic results from the Regge and Coulob regimes respectively in the case of equal 
charges. 

E 



m 




0.0 0.2 0.4 0.6 0.8 1.0 x 1 ' 4 



Figure 6: E/m as a function of J for B\fX/ra 2 — 0.1 (blue squares), B — (red circles) and 
B\fX/m 2 = —0.1 (green diamonds). 

by A J = Bp 2 \^ + . We present our results for the meson spectrum in this case as a function 
of B/m? together with a comparison to the analytic results for the 1 << J « vA and 
J » vA regimes in figure (||) . 

In figure @ we show the energies for B a = ±0.1 as a function of J. In the left panel 
of figure (0) we also show one example of a calculation performed at an intermediate J for 
arbitrary B exhibiting interesting non-linear effects in B. 

In the case of equal opposite charges one has to deal with the additional complication 
that the configuration is no longer symmetric with respect to reflection around the string 
midpoint. In our numerical scheme this means we can no longer assume that the turnaround 



- 10 - 



1.90 



1.85 



1.80 



1.75 



1.70 L 



ByfY 



1.853 
1.852 
1.851 
1.850 
1.849 
1.848 
1.847 



K 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 



0.0000 0.0005 0.00 1 0.00 1 5 0.0020 0.0025 0.0030 0.0035 ,„2 



Figure 7: Left panel: Meson spectrum E/m as a function of B/ni 2 for Jj\f\ = 0.207 for equal 
charges. Right Panel: Meson spectrum E/m as a function of B/m 2 for J /\[\ = 0.255 for equal 
opposite charges. Note that for equal opposite charges the linear Zeeman effect does not contribute 
since the meson has no net magnetic moment, the shift starts at quadratic order in B. 
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Figure 8: String profile r(p) for the case of equal opposite charges with ^ = 0.126 and -M^ = 

point is at p = 0. We therefore must chose both tq and r' Q at p = 0. That is at <jq = r$ we 
impose the boundary conditions 

r(ff ) = r 0j r'(a )=r' . (4.11) 

For every r m one can once more find a B that makes the UV boundary condition true at 
either end of the string. But only for special values of r m will one get the boundary condition 
to be true with the same B at both ends. We therefore use a shooting algorithm where, for 
a given string solution, we step through possible values of the mass to find the special r m 's 
for which the asymmetric boundary conditions are satisfied. We present our results for the 
meson spectrum in this case in the right panel of figure (|7|). A typical string configuration 
for this case of equal opposite charges is displayed in figure @ 
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4.3 J << \/A, equal charges 

An analytic solution to the equations of motion can be found in the case uj — > oo, following a 
procedure similar to the one in [7]. We will see that this case corresponds to J << \/A. To 
simplify the equations, we use the coordinates z = ujz = uj/r,p = ujp and the gauge p = a. 
The equation of motion for z is then 

" A~A*=* (4-12) 



/" - 2 - = (4.13) 



1 + z' 2 z 1-p 2 

We can substitute the ansatz z = ujzm + uz f(p) into (|4.12|) , where zm is the position of 
the D7 brane. Keeping only terms of order 1/ojzdi we get an equation for /, 

pf_ 
1-p 2 

From the boundary conditions for z, we find the boundary conditions /(0) = and /'(0) = 0. 
Then, solving for f, 

/(P) = ^(p 2 + arcsin 2 (p)). (4.14) 
At the boundary p must satisfy 7r p = §j? = —pBu. So, 



VX 1 fl 



P 



2 



, v , ■ -pB (4.15) 

2vr z 2 V 1 + F 2 P V ; 

Substituting the equation for z, we see that the critical value of p at the boundary must 

satisfy, for small B, 

fe- ^ .l.y (4.16) 

Putting all of this into the equations for E and J, including the endpoint contribution to J 
( |4.10| ), using that zm = and keeping terms linear in B, we find that 

£=4^-^ (4.17) 

\/X B , 

J = 4 2 2 " -2 4 - 18 

Az 2 D7 ui 2 oj 2 

Since w is large, J << vA, as we expected. Solving for £7 in terms of J, using ( |2.10| ) and 
again keeping only terms linear in B, 

m A 1 / 4 7T m 2 " 1 ; 

To find the energy of the upper branch, we must integrate from pc to p = 1, and then from 
p = 1 back to p = to take into account the shape of the string configurations as we see in 
figure @). We find the energy of this branch to be, 



m A 1 / 4 vr m 
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So the splitting is AE/m = \/~J\ 1 / i B /m 2 ir, which is consistant with what was obtained in 

(PD- 



4.4 J « vA 5 opposite charges 

For the case of equal and opposite charges, the boundary condition is different for the left 
and right ends of the string. 

tt£| ± = ±pBcu (4.21) 

As we can see from the string profile in figure (Q), on one side of the string, we must again 
integrate from pc to p = 1 to p = 0, as in the upper branch of the equal charge case. The 
other side of the string is integrated as in the lower branch of the equal charges case. Following 
the same method as before, we obtain, 

E = (4.22) 

z 2 m uj 

SO 

E/m = (4.24) 

which does not depend linearly on B. In this case, exact evaluation of the integrals of the 
small J solution without using the additional approximation ( [4.16| ) shows that E does not 
depend on B. To see what this means, consider writing E 2 as a power series in J with 
functional prefactors Ci(B), 

{E/m) 2 = Y,Ci{B)J i . (4.25) 

i 

The result above is simply the statement that C\ is the constant 47T 2 /■s/X. This is in perfect 
agreement with the result of [24] that for a field theory Nambu-Goto flux tube with equal 
opposite charges the energy is completely independent of B. 



5. Meson dissociation in the presence of magnetic fields 

Last but not least, we want to put together the two cases we studied in the last two sections 
and study the strings when both temperature and external magnetic field are turned on. 
There are of course now many different parameters to tune: the temperature, the mass, the 
magnetic field, and the angular momentum. The main point we want to make here is that 
one can map out the properties of the strings as a function of all those control parameters 
using our technique. One particularly interesting quantity we calculate is once more the 
critical angular momentum J cr it(m/T, B /m 2 ) beyond which the mesons dissociate. In figure 
(||) we show J cr it as a function oim/T for B = ±0.157\/Ar 2 . As expected, the magnetic field 
increases the binding energy of at least one state and hence increases J cr u when compared to 
the zero temperature cases. 
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Figure 9: J crlt as a function of m/T for — —0.157 (green diamonds), (blue circles) and 

+0.157 (red squares). 
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